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Abstract. Active galactic nuclei (AGN) are known as sources of high energy 7-rays. The emission probably results 
from non-thermal radiation of relativistic jets belonging to the AGN. Earlier investigations of these processes have 
suggested that neutrinos are among the radiation products of the jets and may be used to discriminate between 
hadrons and leptons as primary particles for the production of the high energy emission. Our calculation of the 
high energy neutrino emission from the jets of AGN is based on a recently published model for 7-ray production 
by a coUimated, relativistic blast wave, in which the spectral evolution of energetic particles is determined by 
the interplay between the particle injection by sweep-up of the interstellar medium, the energy losses through 
radiation, and diffusive escape. It is important to note that the swept-up interstellar particles retain their relative 
velocities with respect to the jet plasma, but get isotropised in the jet rest frame by self-excited Alfvenic turbulence. 
The bulk of the neutrino emission is expected in the energy range between 100 GeV and a few TeV. It is shown 
that the neutrino flux is correlated with the flux of TeV 7-rays. This allows to distinctly search for neutrino 
emissions from the jets of AGN by using the TeV 7-ray light curves to drastically reduce the temporal and spatial 
parameter space. Given the observed TeV photon fluxes from nearby BL Lacs the neutrino signal from AGN may 
be detectable with future neutrino observatories as least as sensitive as IceCube. 

Key words. Galaxies: jets - BL Lacertae objects: general - Gamma rays: theory - Radiation mechanisms: non- 
thermal 



1. Introduction 

To date, more than 60 blazar-type active galactic nuclei 
(AGN) have been detected as emit ters of MeV-GeV 7-rays 
by E GRET (von Montigny et al. |l995l Mukherjee et al. 



1997), and another few BL Lacertae objects have been ob- 



served as sources of TeV 7-rays by ground-based Cerenkov 
telescopes (Catanese & Weekes 1999). In these sources the 
bulk of luminosity is often emitted in the form of 7-rays. 
The 7-ray emission of the typical blazar is highly variable 
on all timescales down to the observational limits of days 
at GeV energies and hours at TeV energies (e.g. Mattox 
et al. 1997, Gaidos et al. 1996). Emission constraints such 



as the compactness limit and the Elliot-Shapiro relation 



(Elliot & Shapiro 1974) are violated in some 7-ray blazars, 
even when these limits are calculated in the Klein-Nishina 
Hmit (Pohl et al. |1995[ Dermer & Gehrels |1995D , which 
implies relativistic bulk motion within the sources. This 
conclusion is further supported by observations of ap- 
parent superluminal motion in many 7-ray blazars (e.g. 
Pohl et al. 11 9951 Barthel et al. |1995|, Piner & Kingham 



1997a, 1997b) 



Most published models for the 7-ray emission of 
blazars are based on inverse Compton scattering of soft 



target photons by highly relativistic electrons in the jets 
of these sources. The target photons m ay co me directly 
from an accretion disk (Dermer et al. |1992|, Dermer &; 



Schlickeiser |1993| ), or may b e rescattered accretions disk 
emission (Sikora et al. 1994), or may be produced in the 
jet itself via synchrotron radiation (e.g. Bloom & Marscher 



1996 and references therein). The criticism of these models 



focuses on the fact that an electron acceleration process 
would have to be very fast to compete efficiently with the 
radiative losses at high electron energies, and it is there- 
fore unlikely that the electrons are the primary particles. 

The radiating electrons can also be secondary parti- 
cles produced in inelastic collisions by primary hadrons, 
in which case neutrinos would be emitted in parallel to 
the 7-rays. Neutrinos can therefore be used to distinguish 
between purely leptonic models and models which involve 
high energy hadrons as primary particles. 

One class of hadronic models is based on proton ac- 
celeration to ultra high energies around 10^" eV, so that 
photomeson production on ambient target photons can 
provide many s econdary electrons and positrons (Kazanas 
Sikora et al. 



& EUison 1986 



1987, Mannheim & Biermann 
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1992 ). Such systems are usually optically thick and a pair 
cascade develops. These models are particularly attractive 
because they allow a possible link be made between escap- 



2 



C. Schuster, M. Pohl and R. Schlickeiser: Neutrinos from active galactic nuclei as a diagnostic tool 



ing ultra high energy hadrons and the cosmic ray spectrum 
at 10^" eV near earth (e.g. Waxman 1995). Steady-state 
calculations indicate that blazars would produce a high 
flux of neutrinos with energies up to 10^^ eV (Mannheim 
19951, Halzen & Zas |1997|, Protheroe |1997D. Note, however. 



that some of these results appear to be in conflict with 
physical upper limits for the total neutrino intensity pro- 
duced by ultra high energy hadrons in AGN (Waxman &; 
Bahcall 1999, Bahcall & Waxman 2001). It is unfortunate 
that time-dependent calculations of proton-induced cas- 
cades have never been performed. We can therefore not 
comment on possible correlations between the neutrino 
flux of a source and its variable GeV-TeV 7-ray emission, 
for the latter is also influenced by the presumably variable 
pair-pair optical depth. 

Two of us have recently published a new model of 
particle energization in jets of AGN, which is motivated 
by the blas t wav e models of gamma-ray bursts (Pohl & 
Schlickeiser 2000 ). There the AGN jet is assumed to be a 
cloud of dense plasma which moves relativistically through 
the interstellar medium of the AGN host galaxy. It was 
shown that swept-up ambient matter is quickly isotropised 
in the blast wave frame by a relativistic two-stream in- 
stability, which provides relativistic particles in the jet 
without invoking any acceleration process. The typical pri- 
mary particle would be a proton with a Lorcntz factor of 
the order of one hundred, distributed isotropically in a 
jet which also moves with a bulk Lorentz factor of the 
order of one hundred. Inelastic proton-proton collisions 
with the background plasma would then lead to the pro- 
duction of TeV 7-rays by Tr'^-decay as well as synchrotron, 
inverse Conipton, and bremsstrahlung emission produced 
by secondary electrons at lower energies. The fully time- 
dependent calculations show that the observed variability 
behavior and 7-ray spectra of blazars can be explained 
with this model. 

In this paper we are concerned with the neutrino emis- 
sion which is produced in parallel to the 7-rays during the 
decay of the charged pions. We shall provide specific pre- 
dictions for the neutrino flux and spectra in relation to the 
TeV 7-ray emission of AGN. These will allow observers to 
a) search for neutrino emission of AGN by using the TeV 
7-ray light curves as indicator of activity, and b) perform 
a decisive test of at least one hadronic model of 7-ray 
production in AGN. 

In Sect. 2 we briefly review the spectral evolution of the 
energetic protons, resulting from the interplay between the 
injection and the losses. We discuss some resulting proton 
spectra and show that the spectral evolution of the proton 
distribution depends on the geometry and the physical 
constitution of the blast wave. In Sect. 3 we discuss the 
various decay modes and present the the source functions 
for the decay processes of the charged pions as well as the 
subsequent decay of charged muons. 

We also investigate the emission of electron anti- 
neutrinos caused by /3-decay of neutrons produced in 
p+p p + n + X reactions. We show that the power radi- 
ated in thus produced emission is about two orders of 



magnitude smaller than the power radiated in the form of 
electron (anti-)neutrinos arising from the decay of muons. 
In Sect. 4 we discuss the spectral evolution of the neu- 
trino emission, induced by the proton spectra presented 
in Sect. 2. We show that (muon)-neutrino emission at en- 
ergies around 1 TeV is correlated with TeV 7-emission. 

2. The energy distribution of protons in an AGN 
jet 



The AGN model of Pohl & Schlickeiser (2000) describes 



the jet as a plasma could consisting of electrons and pro- 
tons, that is assumed to have a cylindrical shape with 
thickness d, where d is small compared to the radius R 
(see also Fig.0). The jet is moving with a bulk Lorentz 
factor r through the ambient medium. Viewed in the rest 
frame of the jet plasma cloud, the interstellar medium 
forms an electron-proton beam of density]^ Ui — T n* , 
which is small compared with the plasma density rib of the 
jet cloud. Using kinetic theory Pohl & Schlickeiser (2000) 
have shown, that the incoming particles of the interstel- 
lar medium quickly get isotropised by self-excited Alfvenic 
turbulence, but retain their relative velocities with respect 
to the jet plasma. The jet plasma thus becomes enriched 
with relativistic particles, of which only the protons are of 
interest for us, for they can produce neutrinos. 



n 



r 




d 



Fig. 1. Sketch of the basic geometry. The thickness of the 
channeled blast wave d, measured in its rest frame, is much 
smaller than its halfdiameter. The blast wave moves with 
a bulk Lorentz factor F through ambient matter of density 



^ All quantities indexed with 
(galaxy) frame. 



are in the laboratory 
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Fig. 2. A typical spectral evolution for protons in the blast wave, responsible for the neutrino emission. In this 
example the following parameters have been used: The radius of the plasma disk is i? = 10^'' cm, the thickness of the 
disk amounts to d = 3 • 10^'^ cm and the initial Lorentz factor is Fq = 300. The constant densities inside and outside 
the jet are rib = 5 • 10^ cm^"^ and n* = 0.2 cm~^, respectively. Note that the mentioned time refers to the observer 
frame and therefore depends on the viewing angle ^, which we choose to be 0.1° in this example. 



This has two consequences: 
a) the system is not stationary. For reasons of momentum 
conservation the Lorenz factor F decreases at a rate 



• ^ TT n* c nip 2 _ -|n3/2 
M{t) ^ ' ■ 

Note that the relativistic mass 



(1) 



M{t)^ jmpN{^,t)d^ (2) 

is obtained from the particle spectrum in the jet plasma 
N{-f,t) = V{n{-f,t)+nbS{-f-l-e)). (3) 
in which, in contrast to the nomenclature used in Pohl and 



Schlickeiser (2000), we have have included the background 



kT 



< 1. 



plasma with a formal thermal energy e 
b) relativistic protons are injected in the jet plasma at a 
rate 



5(7, t) = TTR'n*c^/TW^5{j ~ F(t)) (4) 

where F(t) is determined by Eq. ^ The temporal evolu- 
tion of the proton spectra N{'y,t) can be described by the 
continuity equation: 

dNij,t) , d{jNh,t)) , iV(7,<) , 7V(7,<) 



dt 



97 



T. 



N 



Sil,t) (5) 



We consider continuous energy losses from clastic and in- 
elastic scattering at a rate of 



7 



3 • IQ-^^nt 



7 



7- 10" 



-16 



1 



rib- 



( 7-1)^ 
7 + 1 



(6) 



Additionally, one has to consider catastrophic particle 
losses arising from diffusive escape at the timescale of 



T. = lt c.l.6710--'^'r< 



A/3c 



(7) 



where A is the scattering length and /3 — \/l ~ "f ^. 

Neutrons produced in p n reactions are likely to 



10 



15 , 



s for 



escape the jet at a timescale of Tjv — 3 
7 > 1 (cf. Pohl and Schlickeiser |2000D . 

Neglecting energy losses and additional particle losses, 
the continuity equation (|^) becomes solvable and we get a 
(7^ — 1)~^/^ spectrum for the particle distribution function 
Af(7,i). However, more realistic cases require the use of 
numerical methods. 

We use Eq. (|l|) to obtain the time r* at which the blast 
wave would be observed with a particular Lorentzfaktor 
F. 



M(r') 



c(r'2- i)(3/2) 



r'(i-/3cosr) 



(8) 



where C = tt n* c nip and f3 — \f\--V^ . For the pur- 
pose of transformation to an observer frame one should 
bear in mind that the energy scales as well with the factor 
r(l-/3cos6l*)-i. 

In Fig. H and ^ we calculate typical examples of result- 
ing proton distributions -/V(7,t), describing the number 
density of protons in terms of the Lorentz factor 7 w.r.t. 
the blast frame. Starting with an initial Lorentz factor of 
Fq = 300 we show the distribution of protons for several 
instances of observed time corresponding to the reduced 
bulk Lorentz factors of the blast wave. According to Eq. 
|l| the Lorentz factor is monotonically decreasing. In all 
cases the sweep-up occurs at the present Lorentz factor in 
accordance with Eq. ^ 

As mentioned before, the overall behaviour of the par- 
ticle spectra is mainly determined by the competing modes 
of sweep-up and cooling of the particles. In Fig. || the 
swept-up particles cool down faster than the blast wave 
decelerates, whereas in Fig. ^, down to a Lorentz factor 
of F = 85 the cooling is slow compared with the deceler- 
ation of the blast wave. Afterwards the deceleration rate 
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Fig. 3. The proton spectra calculated with another set of parameters. Here we use R = 2 ■ 10^^ cm for the radius 
of the plasma disk and d = 10^'' cm for the thickness of the disk. The initial Lorentz factor stays at Fq = 300. The 
density in the jet has changed to Ub — 10^ cm^'^ and for the density outside the blast wave n* = 1.5 cm~'^ is assumed. 
Here the specified time refers to an viewing angle oi 9 — 2°. In the bottom row the same spectra as in the top row are 
shown but with higher resolution. 



of the blast wave is overtaken by the cooling rate of the 
particles. Then the competing mode of particle cooling is 
the more effective process. Determined by the initial pa- 
rameters d, R, rii* and rif, the spectra show a very different 
behavior. The time specified in the pictures refers to the 
observer frame and therefore depends on the viewing angle 



3.1. Pion decay 

The pion decay muon neutrino source function {Ei, , t) 
describing the energy spectrum of the neutrinos produced 
in the jet via pion decay can be written as 



(9) 



3. The production of high energetic neutrinos 

The energetic protons in the jet plasma produce pions 
in inelastic collisions with particles in the background 
plasma. For the calculation of the pion source spectra we 
use the Monte -Carlo model DTUNUC (V2.2) (M ohring 
& Ranft 119911, Ranft et al. 11994, Ferrari et al. |l996a|, 



Engel et al. 1997), which is based on a dual parton model 



(Capella et al. 1994). Here we are concerned with high en- 
ergetic neutrinos resulting from the decays of pions tt^ — > 
11^+ Vfj.i'i'fi) and subsequently /i^ — > e^-|- z^e(i'e)- 



where is the pion source function and 

is the normalized neutrino spectrum describing the two- 



body decay in t he res t frame of the plasma (Stecker |1971 
Marscher et al. |l980| ). Here E° ~ 29.79 MeV denotes the 



energy of the emitted neutrino in the CMS. The lower 
limit of integration is given by 



2E„ 



2E?/ 



(11) 



The pion source function giving the amount of pions per 
cm'^, GeV and s can be calculated as an integral over the 
spectrum of energetic protons and the differential pion 



C. Schuster, M. Pohl and R. Schlickeiser: Neutrinos from active galactic nuclei as a diagnostic tool 



5 




Fig. 4. The production rate of muon neutrinos resulting from the various decay modes calculated for one proton of 
kinetic energy 10 GeV and 100 GeV, respectively. The emission rate does not change in the energy regime of a proton 
up to a range of 100 GeV. Additionally we show the emission of 7-rays resulting from the decay 7r° 27 displayed 
with the solid line. The emission from the positively charged pions and muons is slightly higher than the emission 
from the negative ones. The respective pion decay is more effective. 




log(E [GeV]) log(E [GeV]) 

Fig. 5. For a proton with the same kinetic energy as above the electron neutrino emission is shown here. The production 
ofVe and Ve resulting from the muon decay is displayed by the solid line. The energy of anti-electron neutrinos produced 
by /3-decay is about 2 orders of magnitude smaller. The complete rates of muon neutrinos and 7-rays are also depicted 
for reference. 
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(1 + ^T+|;)3(4C- 3(1 + /3J)C') 



production cross section which we derive using the 
Monte-Carlo code DTUNUC. 



(12) 



for 



/3p - . - 1 



Here p = Vp/c denotes velocity of the protons. It is con- 
venient to perform the integration over the pion Lorentz -^^j^gj-g a 
factor in Eq. first. Hence 



< c < 



(18) 



and /3p = ■u^/c. Here J7vhfi,x) 



denotes the Lorentz factor of the parent pion: 



— — —d-/Tr d-jp (13) 



Si' 



(19) 



is the differential production rate of particles per cm^, 
GeV and s. 

3.2. The neutrino production via decay of charged 
muons 

In the decay of the pions charged muons emerge producing 
neutrinos in a second decay /i^ — > + V^{v^) + VeiVe)- 
The muon decay source functions for muon and electron 
neutrinos Q'^{Ei,,t) are: 



oo i 



dxdlt. 



(14) 



7m. » 



with the neutrino spectrum f{E^,^^) in the rest frame 
of the plasma describing the three-body muon de- 
cay (Marscher et al. 1980, Zatscpin & Kuz'min 19621 ): 



RE.) 



16 7m'( — 
7m 



-(3 + /3/)C)C' 



for < C < 
and 



1-/3^ 



(15) 



f{E.) 



(1 + 



(1+/3m) 



r(^-3(l + /3^))C^) 



for 



1-/3m 



< c < 



111 c2 



(16) 



for muon neutrinos and electron neutrinos Ve- The spec- 
tra for the production of anti-muon neutrinos and anti- 
electron neutrinos are calculated by: 



/(£;,) - 32 (_ 
7m 



2(3 + /3,,')C)C' 



1 



for < C < 
and 



1-/3m 



(17) 



where ei = 3.68, £2 = 3.55 7^ and x is the CMS interac- 
tion angle cosine. 

The lower limits of integration of dll) are: 



E, 



and 

Xrriin 



(20) 



£1 > £2 



The calculation of the source function (57r(77r(7MJ x): ^) is 
based on the Monte Carlo code DTUNUC again and is 
calculated analogous to the previous case with formula 

3.3. Neutron (3-decay 

The emission of electron anti-neutrinos resulting from the 
/9-decay of secondary neutrons n ^ p + e + Ve is described 
by the electron /3-decay so urce f unction QiJe (E^, , t) calcu- 
lated analogous to Jones (1963) treating the neutrino as 
a particle without mass: 



1 ffiE^u) f Qniln,t) 



E" 



1 



: d7n dEl 



(21) 



with the upper limit of integration Q — m„ (? — nip (? 
and the lower limit of integration 7mm given as: 

1 E E^ 

Imin = 7: (-^ + -J=r) (22) 



The CMS normalized antineutrino spectrum is given as 
(Leon |1973| ): 



f{El) ^C{Q~ ElW[Q-E^,Y-mlc^ Ef 



(23) 



Again E'^ denotes the energy of the emitted neutrino in 
the CMS. C is the constant for normalization and is 1.2 
MeV~^. For the description of the neutron source function 
Qnilmi) we use a discrete sum of intensities of nionocn- 
ergetic neutrons, each one modeled as a ^-distribution. So 
one assumes that the thermal protons take part in the 



C. Schuster, M. Pohl and R. Schlickeiser: Neutrinos from active galactic nuclei as a diagnostic tool 



7 




S- 

O" 

H 



-10 12 3 4 
log(E [GeV]) 




H 



-10 12 3 
log(E [GeV]) 






-9 




a 








-10 




H_ 








-11 












-12 






J I I L 



J L 



log(E peV]) 



> 

a; 
O 




> 

O 



log(E peV]; 




log(E peV]; 



Fig. 6. The evolution of the muon neutrino emission resulting from the proton spectra of Fig. S. The viewing angle 
is 6* = 0.1° and the redshift of the AGN is 2; = 0.5. All the other parameters are chosen as in Fig. g. In the top row we 
see the F^, spectra and in the bottom row the v F^, spectra. The spectral evolution of the 7-ray production spectra is 
also shown for reference. 



collision process as well as the relativistic ones. The ap- 
proximation for the source function Qn{ln,t) is: 



[5{ln - ^) + 5(7„ - 1.1)) d7p (24) 



with T = 3 • 10"'^^ ''h'^ ■ second ^- function with the 
argument 1.1 accounts for the thermal protons. Again we 
change the order of integration: 



00 Q oc 



7„ 



d'jn dE° d'jp 



(25) 



with _D(7„) — 5(7,1 — ^) + (5(7„ — 1.1). To obtain the results 
in the observer frame one has to consider the transforma- 
tion: 



q{e4e;),t*) 

r2 (1 - /3cos6i*)2' 



where all quantities indexed by 
observer system. 



(26) 



are calculated in the 



4. The evolution of the neutrino emission spectra 
in comparison with 7-ray production 

In this paper we focus on the high energy muon neutrino 
emission resulting from the decay mode of pions tt^ — s- 
fJ.^+ v^i,{p^) and subsequently /i^ +v ^{v^)+ Ve{pe)- 

In Fig. H and Fig. |^ we show two examples of spectral 
evolution of total muon neutrino emission calculated with 
the model proton spectra seen in FigH and in Fig.^. We 
assume a constant density of the background plasma. The 
production rate of 7-ray resulting from the decay 7r° 27 
is depicted as well for reference. The bulk of muon neu- 
trino emission lies between 100 GeV and 1 TeV and fol- 
lows strictly the evolution of 7-ray production. All spectra 
are calculated using the Monte-Carlo code DTUNUC de- 
scribed earlier in this paper. Emitted 7-ray photons may 
be absorbed by collisions with infrared background pho- 
tons in the intergalactic medium, but neutrino emission is 
unaffected by any interaction. So for known 7-ray produc- 
tion the neutrino emission may even exceed the calculated 
rates. 

Despite large production rates the detection of high 
energy neutrino emission requires large water or ice tele- 
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Fig. 7. The second set of parameters used in Fig. || results in these spectral evolution. In this example we use the 
viewing angle 9 = 2°. The strong dependence of the observed emission from the angle 9 is due to the high Lorentz 
factors in our model. This is responsible for the strong rise after 1000 h. The set of parameters is the same as in Fig. 
S and the redshift is again z = 0.5. 



scopes, based on detection of Cerenkov light (for example 
the AMANDA experiment at the south pole). In this sec- 
tion we want to estimate whether or not neutrino emis- 
sion from blazars will be observable, given their TeV 7-ray 
fluxes. 

Suppose a TeV source is observed with a 7-ray flux of 



~ 10"^° ph. cm 



at a few TeV 



(27) 



As we have seen, this implies a similar flux of muon neu- 
trinos at a bout a TeV. Following the treatment by Gaisser 
& Stanev (1984), we can calculate the differential flux of 
muons produced in deep-inelastic interactions of the in- 
coming neutrinos with the ice material as 



dE„ 



10 



-19 



a 



GeV^^cm-^sec-i (28) 



when the deep-inelastic cross section a, the muon energy 
loss rate a, the Avogadro number iVAv, and the neutrino 
flux according to the 7-ray flux (|2^) are inserted. The 
muon spectrum is flat for i?^ <C Ei, and the mean energy 



1/3 ~ 0.3 TeV. Then the total muon flux is 



E,, 



approximately 



— dNa 

' E, 



dE„ 



3 • 10 



- 1 7 ~9 —1 

^ ' cm ^ sec 



(29) 



The number of neutrino-caused event now depends on the 
effective area of the detector A^s, and the observing time 



tobs- The effective area of the AMANDA detector, which is 
operational to date, is about AcSj^ ~ 10^ m^, whereas the 
projected IceCube experiment may have A^Si — 10^ m^. 

The number of events from a TeV 7-ray source with 
photon flux (|27|) in the two detectors is then 



iVA ^ 0.01 



I Month / 



\ Month / 



(30) 



Equation ( p(]| ) makes clear that the observation of neu- 
trinos from AGN is, if at all, only possible with future 
neutrinos observatories like IceCube. 

Is it possible to identify an AGN neutrino signal in the 
data of IceCube? The actual detection of some neutrino 
events in the detector is not sufficient for them being iden- 
tified as coming from an AGN. We need to observe a local 
excess of events on top of the nearly isotropical distribu- 
tion of events caused by atmospheric neutrinos. The inten- 
sity of atmospheric neutr inos a t 1 TeV is about (Volkova, 
Zatsepin and Kuzmichev 198C) 
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-7 —2 —1 —1 

cm sec sr 



at 1 TeV 



(31) 



It varies by a factor of five between vertically and horizon- 
tally incident particles, but can be taken as locally flat. At 
muon energies around 300 GeV the muon follows the flight 
direction to within slightly more than one degree. No mat- 
ter how well the detector can reconstruct the muon flight 
direction, the uncertainty in the reconstructed neutrino 
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Required event rate 




10 12 



Fig. 8. The required rate of AGN neutrino events in 
IceCube for a 3a excess per resolution element as a func- 
tion of the observing time. The neutrino rate expected for 
a TeV 7-ray source with a 7-ray flux as in Eq. ^ would 
be one per month. Therefore an observing time of nearly 
one year would be needed to produce a 3cr signal. 



flight direction will in any case be worse than one degree. 
For our purpose it may suffice to assume that the neu- 
trino incidence direction can be reconstructed to within 
two degrees, which leaves some room for uncertainties in 
the muon event reconstruction. We can then define a reso- 
lution element, fl ~ 3.8 • 10""^ sr, and the number of events 
of atmospheric neutrinos per resolution element A'atm for 
Icecube as 

Natm = ^ ^cff: t ~ 4 at 1 TcV (32) 



which is four times the expected event number of neutrinos 
from an AGN with 7-ray flux . 

The expected number A^agn of atmospheric neutrino 
events per resolution element would be well defined, for 
we can use thousands of resolution elements to determine 
it. Given an observed number of events, A^obs, we would 
deduce the number of events due to AGN neutrinos as 



A^AGN = Aob. - Aa 



A, 



obs 



_4 Job^\ 
V Month J 



(33) 



In the null hypothesis A'obs follows a Poissonian distri- 
bution. 



PiNobs,Natm) 



^ atm 

N J, I 

ohs- 



exp{-Natrn) 



(34) 



Equation (^3|) allows us to determine the event rate of 
AGN neutrinos required to provide a signal of specified 



significance by appropriate summation of the Poissonian 
probability distribution ( p4|). As an example we give the 
result for a 3cr excess in Figpl Please note that the signifi- 
cance applies for an excess in an predetermined resolution 
element. 

Apparently much more than 10 events are necessary to 
detect signals from AGN. By approximating the width of 
the Poissonian distribution by the square root of its mean 
we can therefore derive a simple formula to calculate the 
signal to noise ratio as a function of the exposure time i, 
effective area A and the 7 source flux A^^ from AGN: 



VlO^m^y V Month/ 



1/2 



(35) 



Mkn501 has been observed to emit a 7-ray flux of 
0.5 • lO^^'^cm"^ s~^ at a few TeV for about half a year 



(Quinn et al. 1999 ). Mkn421 has been similarily active in 
the first haff of 2001, as discussed at the ICRC 2001. So 
this indicates the fiare states of many AGN would have 
to be co-added to expect a meaningful neutrino signal in 
the IceCube detector. Nevertheless, because forthcoming 
TcV Cherenkov telescopes will provide much better sam- 
pled light curves of many more AGN, the detection of a 
neutrino signal from AGN is still possible. 



5. Discussion 

We have calculated the neutrino emission resulting from 
jets of AGN, under the assumption of a channeled blast 
wave model. Therefore we have investigated the decay 
modes of pion and subsequent muon decay. Neutrino emis- 
sion resulting from the decay of secondary neutrons has 
been described as well. We have shown that neutrinos re- 
sulting from /?-decay of neutrons possess an energy about 
two orders of magnitude lower than that of the other decay 
channels we considered. The bulk of the neutrino emission 
is expected to be in the energy range between 100 GeV 
and 1 TeV. 

The emission rate resulting from a proton pair in the 
energy regime between 1 GeV and 100 GeV does not vary 
much. Therefore the time dependence of the emission spec- 
tra is completely determined by the changing spectra of 
incoming protons. We additionally discuss some spectra 
resulting from the channeled blast wave model of Pohl 



and Schlickeiser (2000) and observe that their shape is 
determined by the ratio between the cooling rate of the 
particles in the blast wave and the deceleration rate of the 
blast wave itself. 

Neutrino emission should be correlated with the emis- 
sion of 7-rays. It allows to distinctly look for neutrino 
emission from the jets of AGN by using the TeV 7-ray 
light curves to drastically reduce the temporal and spa- 
tial parameter space in the search for neutrino outbursts. 
Given the observed TeV photon fluxes from nearby BL 
Lacs the neutrino flux can be detectable with future neu- 
trino observatories. 
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